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Contributions to the theory of diffraction of electromagnetic 
waves by spherical particles. II 


By T. LiunccrEN 


With 2 figures in the text 


1. Introduction 


In a previous paper! a new method of attacking the numerical evaluation 
problem of diffraction of electromagnetic waves by spherical particles has been 
presented by the author. In the following it will be referred to as A. Here, 
the phenomenon in the forward direction will be studied in some detail. In A 
only the term interpreted as describing the diffracted radiation was computed 
for that direction; the result (eq. (A 102)): 


wt . —tkor J 
Es,, =1 0” cos ee a 
Egy, = t8 PH; | 
y=090 


agreed with results won by elementary methods. 

The practical consequence of (4 102) is a much used method of determining 
the size of spherical particles, e.g. in mists and clouds, by observations on the 
minima of intensity in the forward direction.2, These minima lie in directions 
determined by the zeros of the Bessel function®: 


J1(99) =0; (1) 
or explicitly, for the first few: 
Yi = 0 9, = 3.832; yo = 0 Og = 7.016; y3 = oe O3 = 10.173. (1 a) 


For large particles this theory is quite satisfactory, but both Mrckxr and 
Kou er (loc. cit.) emphasize that for decreasing radius of the sphere it must 


1 LagunGGREN, Ark. Mat. Astr. Fys., 836A, n:o 14, 1948. 

* See e.g. K6uLER, Untersuchungen tuber die Elemente des Nebels und der Wolken, Diss., 
1925; Mrcxr, Ann. d. Ph., 61, p. 471 & 62, p. 623, 1920. 

3 In elementary theory we have (cf Handbuch der Physik, XX, p. 179) sin © instead of 
©. As is pointed out in A, pp. 33 & 34, the argument @ © is a consequence of a better 
approximation of the Legendre functions. 
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meet with grave objections. Mecke’s theoretical investigation, using geometrical 
optical methods, does not give an ultimate answer to the question when eq. (1) 
can be trusted. Especially, Kéhler uses (1) far below the vague limit given 
by Mecke (@ ~ 50—70; our notations), and obviously this causes no difficulties. 
In this paper we will reexamine the problem, starting from an exact solution 
of Maxwell’s equations. 

As has been mentioned in A, this investigation implies a discussion of the 
reflected and the refracted radiation, that is of the terms Cy, and C, for 
Q@~0; the application of the saddle point method encounters the difficulty 
with a saddle point in the vicinity of the limits of integration. These two 
cases are, however, not equivalent; the treatment of the reflected radiation (C3) 
will be more intricate, as we will then receive the essential contribution from 
the region of the integration interval in which the asymptotic behaviour of 
the cylinder functions changes. 


2. The refracted radiation 


We start with (A 25), applying (A 20) to get a description valid for 0 = 0: 


es “(4-220 N 1 
Ci = nye fe Fea Ty las (ne) — nase 


- [2 (o) C2’ (no)]-2 PY” (cos O) sm OdE. (2) 
The term C3) is obtained through the exchange: 
[a3 (ne) — nage: (0)]-? > — [nqiz (me) — ge:(0)]*; 


and the terms with an upper index (1) through another exchange: 


a, P’? (cos @) 
PS? > = 3 
?” (cos 9) sin O ers 
Hq:s (49) then yield the field components: 
vr ___ U COS ¢ Ae —tikor (QV __ yor) 
a ea | 
v0 sin Wee 
W Degen — e Veo? (OY? 
/ kor (C4 CR > | 


The Debye approximations of the cylinder functions (cf. eq:s (A 26)—(4 30)) 
and recurrence formulae for spherical harmonies (ef. (4 32)) now change (2) into: 


e sie cos? O = 
COP 1) iy} 2 ‘XN A N 
yal ee O41 + ane oa | 1) | 


sin” T, SiN T. cos Ty 
(n sin tT, + sin T.)” 


. 4 2%0 [sin 7—n sin 72+¢ 3—-%]1—-AaN z ¢ 
Su:(O)e i pie Se Sri bl hr Races OS | 
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where: 
a 
Su:(0) = cos © Ps (cos 9) + P:+1 (cos al;| 


(3) 
| 
COS Ty = N COS Tz = ——*- 


@ 


bole 


For Su:(@) we use an approximate integral representation, easily deduced 


from (A 97): 


*rola 


2 
sus(6)=*(1+ 2] | cos (0 O cos T, cos yp) dy; (4) 


0 


valid for small values of 9 ( < 4 and especially appropriate for small values 


> 


of €. This is just what we want as the saddle point is situated near tT, = 


i.e. €=0. We further put: 


no] a 


ne ates oll = Slee N 
2 cos O 5 
a 


6 
1% at most. The introduction of these approximations transforms (2a) to: 


which is also a permissible approximation for @ < —; the error is then about 


a4 
2 2 
2 ; sin? T, SIN T. COS T 
= N 1 2 : 
hee fay | arg 
% 2 (nm sin Ty + SiN T9) 
0 


N =— 
0 
[ere SP + ef? OP] dt,; (5) 


Lie st 2 
on + 0% = "2 (1+ 5) 
ait 3 


where: 
G.”_ = cos t,[— 22 N — 2%, + 27, + O cos y|-+ 2 sin t, — 2” sin 7%; (6) 


This is essentially the same exponent as F‘”_ occurring in the calculations of 
‘t1 for 0~0 or x (see A. p. 27 ff.), and of course the calculations given there 
can be reused. We have a saddle point if 


0 
%—%=tZesyp—aN; 


For N =O we get an acceptable value for t,, close to =, whichever sign we 
choose before @. With the plus sign (the term (0; +)) we get a value just 
above ae in the term (0; —) we get a value below = Thus we must take 


2 
into account a saddle point situated just outside the interval of integration. 
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From (477) and (478) applied to the actual case, we get for the saddle point 
T,, in the term (0; +) and for t,, in the term (0; —): 


10 


asin" * 


— C08 T,, = 608 T,, = = aaa = (7) 
Zz: COS > wat 1 


Remembering that O is a small quantity we can write this: 


Pg ee Bes y (7’) 
2 p 


The denominator » means: 
p=|/ SBP Te Seah e APARS gman ate (8) 


In (7) we have excluded the case n ~ 1. We now expand the slowly varying 
part of the integrand and the exponents G”_ in series from the saddle point. 
With the temporary notations: 


Uae Tel 


oh -. (9) 
Usa Gham CAs J 
the result is written 
sin” 1.810 To COST _ 1 Sy OSGO8 Yo). 4 uaewen chan 
(n sin tT, + sin T)? (nm ne is 2p sp. = “a 
) 
1 n@ cos p 
Sa MRTG cs Ce bees fe 
(n + 1)? 2p 
and 
Ge es) = — 2 (a1) - ga = ) @ cos? , ah ee | 
4 p* 
2(n—1) (n—1)(1 + 2) : “ 
2(n — n— n—n 
(Fe e2\ ve Ms sz Pua 2 es ae aes 
(Tio) Ps 3 ne 0" cos” p | 
for both s=+ and s=—. 

We now perform the inner integrations in (5), omitting all terms but (0; +) 
and (0; —). The appropriate path of integration consists, for each term, of a 
straight line of the type (471): 

Re (t — t19) = 1 m(t — Fp); (11) 
from tT — Ty =— (1 + 2) co down to the saddle point, and of the part of the 


Lh 


real axis between the saddle point and t,= 5 - The contributions from the 


path (11) are obtained from the formula: 
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0 et m+1 
vate — 1)” ath ear 1 
| y™ eer dy = oar (=) : r("3 Sem te eoe 0). (12) 
—(1+1) 00 : 3 z 


Denoting by R the contributions from the remaining piece of the integration 
path we get: 


a 8n2 io e2ig@—Y Ee? 
Oe GSP eee See og EE ee = lh ' 
AltOC a1 a(n + 1) (1 = 3 
exp ( Zee ) 9p cos” y 
dy +h. (13) 
Shae a 
J om y_ eo ) 6 cos? y 
Pp 


In the denominator of the integrand the second term can be neglected. The 
zee : a oe 8 
error arising from this exclusion is less than 2 % for @ > 45 and O< a In 


the following we will use the notation: 


pa eee. (14) 
and thus the expression for C4 + C41 is abridged to: 
2 2 2 
Bisciat ig-esee 8 (ie 8) fetoorvaysn as) 
0 


4 : 8 : 
We note that £ is no large quantity. In particular, for 0 or n> 4/3; 
o > 20, we get: 


B<16; (14’) 
For R we have: 
4no? e2iem—)) =| [ iy ee 
i 3 ? 1 eho 27 3 cos? y. d 
ee EE Tp PM i iy . 
: b cos b cos wy 


E b cos p | eeu’ dy —2 f yeien ay| 
0 0 
The new notations mean: 
ay — 4G” (To) 


nO 
sre 
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We have used the fact that the terms (0; +) and (0; —) give the same con- 
tribution to R. 

The second of the inner integrals can be calculated at once; this part of R 
becomes: 


ne 


: 
Batter een “| / ae cae 
ae cae. 1 iB cosy — 1) dy, 
moe ae aye oe 
0 


Denoting the rest of R by R’ we can now write for (13’): 


(2) a) Zio (n—1) n* 1 oe" R’ 15 
e) — 97 ne —2i0 (n—1). ; aes SF tee : 
Gait Ca 2aee (n —1)(n + 1)? ( a 7) 
where: 
n G6 COs w 
22 eS 
ae 4n° 0° O sete). (1 ag =) fer Pes ce eo 2i 8 cos? yr [ efem dy, (16) 
a(n + 1)? ; 
0 
The inner integral is essentially a Fresnel integral: 
nO COS Ww 
20 V2B-cos w 
a ea i ae 
| etal dy = —— [ewan 
f Vom. 
0 
on which we use the formula due to LomMmst!: 
Voz 
| edu =V2z2e8(2); (17) 


0 
where: 


Ve 


@= 3-19 Pim 8-1 |/ Homes 7} 


is a function whose values are readily obtained from the existing rather ex- 
tensive tables of spherical Bessel functions. For values of z considered in this 
paper the series in (17’) converge rapidly. The introduction of S(z) takes R’ into: 


1 


2 
) 2 3 2 
aro “¢ C a Ma —2io(m—-1),. 2 2" —i8 cos? y 
R 2 (1+ ieee ee P é S(B cos? py) dy. (18) 


* See e.g. Watson, Theory of Bessel Functions, p. 546. 
* Tables of Spherical Bessel Functions, I & II, M.T.P., 1947. 
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We can at once see that for very small values of @ this integral tends to 
Xt 


4(n — 1)? 
To get a description for a wider range of @ we proceed as follows. Numer- 


ical checking confirms that for 9 @ <8 and @ > 20, we can use the approxima- 
tion: 


[Js (z) — tJ s (z) — Js (2)]. (19) 


For 00 <4 we may even drop the last term. We denote the contribution to 
the integral in (18) from the first two terms by J, and from the third by Jp. 
We get: 


Zz, == oe | |e COs? yr et cos’ wy. a (6 cos” p) YP | 
re peosty |°¥% | 


2 
I,=—= | ar oes [sin cos? y) + 
2 pp / (6 y) 


(20) 


3 cos (8 cos? y) 3: sin (8 cos? p) 
2 2 dy 
B cos" p fp" cost p 


as p?P is independent of y. In order to be able to compute these integrals 
without resorting to numerical integration we now suppose f to be independent 
of y. This must be considered tolerable as the dependance on wy is weak, fh 
lying between the limits: 


—1)02 C) 
cae a) 
8 (nt —2n cos —> + 1) 
We will use the value for y=0: 
n(n 1) 6? 
pr 22" e 21) 
8 (n?—2n cos 5 + i) 
It is now convenient to use the expansions: 
y= Dan 008 (m + 3) 2; | 
" (22) 
cos eee Pate oe: Hast 
Z Z pees | 
where: 
y SHY eee ae | 
bm rs Si(m + 3)a—Si(m—})x (22’) 
m + 4 t 
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They transform J, and J, into series of integrals of the type: 


he — iWon CO a (23) 
0 


which is easily brought back to integral representations for the Bessel function 
of the order zero. We get: 


TOs ¢ 
I, = 5 E*Io(2). (24) 


Applying this result, introducing the new expressions for J, and J, in R’, and 
for R’ in (15) we obtain: 

(2) <i) * ‘i Ge 

OF; + OF = 200 e 2tem-)). (1 } 3 ) 


where: 


_g2mt8, _ (2m+3 
v Py eee 
+e 4 o( 4 pb | 25’) 


ee .2m-1 , TIS need 
ta = @*8 T(B) —1—3 D bn ba “Io (7 —*e) = 


m =0 


gp 2mt3 5 2m+3 
etre 


Performing the same calculations for the other C, terms we see that for 
O~0 and with the accuracy of (25): 


Cee, On ee ae : 
Cc. = os pe 152. (26) 


Therefore, (25) at once yields the expressions required to give the field com- 
ponents : 


a Ey nv oe? 
E54 = — tos e7 oA Sy ee was : | 
ene ap 2 (n — 1) (nm + ip | eye (27) 
For very small angles @ we get the considerable simplification : 
En = Eo ios pe ltteama ic Seats oe OP (2) 
ia (n — 1) (n + 1)? 4(n— 1)® : 
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We can check our expressions (27) by comparing them to the formula given 
by Mecxe.t For 0 = 0 we get an intensity: 


— = i ns . 2)’ 28 
to (n — 1)? (n + 1)4 (" : (28) 
whereas Mecke gives” 
4n” a\? 
1}0, Mecke = (n = 1) (n + 1)4 ; (“) ; (29) 


The comparison between the two expressions must favour (28). The original 
expressions for the C, terms (cf. eq. (4 25)) show that the field components 


1 ; é t 
tend to zero as 7 for n — co; this accords with (28) but not with (29). The 


error in Mecke’s deduction is, in fact, easily detected. Setting up eq. (17) p 
631, he uses the Fresnel formulae as though the two refractions were identical. 
They are not. Changing Mecke’s expressions accordingly? we get, in eq. (18) 


4n 
Pp. 631, a factor poi? instead of ———-,; the additional n in the numerator 


4 
+ 1) (n + n +1) 
then transforms (29) to (28). 

For 0 +0 other differences between our expressions and Mecke’s seem to 
occur. From his eq. (21) on p. 632 there follows, for small values of @, an 
intensity dependance on » and O only, whereas our formula seems to give a 
dependance on o as well. A closer investigation shows, however, that there is 
no difference of this kind in the two representations. A numerical computation 
of x tells us that the intensity is almost constant for small values of @; so 
does Mecke’s formula. In fact, for small values of 0, it is easily proved that 
(27') might be given as a result of Mecke’s expressions as well (if, of course, 
the error pointed out is corrected). 


3. The reflected radiation 


The additional complication already mentioned and arising in the treatment 
of the terms Cy, must be considered fatal to the method used until now. It 
-eannot but meet a failure when the asymptotic expressions are unreliable in 
the vicinity of the presumptive saddle point. Our defense in A, pp. 8—9, of 
the application of such expressions in a »wrong region) does not essentially 
change this situation. What we said there can be upheld only if the change 
of the phase for |] ~ is rapid enough to cancel the major part of these terms; 
this must be the case both in the exact series and in our in this case quasi- 
asymptotic substitute for it. 


1 Mecxs, Ann. d. Ph., 62, pp. 631 & 632, 1920. 
2 We have corrected an obvious inadvertancy of Mecke in passing from eq. (21) to (22). 
$ Eq. (17) will then run (Mecke’s notations) : 


4n (» cos ‘ -1) (n— cos ®) 


a,—a" (n?—1)? - 
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Realising this breakdown of the method the author has nothing better to 
suggest than a direct summation. As the available tables of the functions in- 
volved do not allow us, for large values of 0+, to begin with the rigorous 
expressions, we must iuiroddee asymptotic formulae in the same way as in A. 
We start with the form (418) for the Co, terms and apply (A 20): 


<2) a I iz $ oe fh P; Dl O; 
Cro = ids xo2' Pi” (cos @) sin ao 
X=A or B. 


The terms with an upper index (1) are obtained through the same exchange 
as in (2). The new notations mean: 


Peo = Hite) Qin ' 

: Ais (e) [er (e) ? [git(n 0) — naei(o)]’ (31) 
F _ Hits (o) | Sa a ia: 

Be =~ Fs (0) [ (@Elnau(me) — az(ol 


We now introduce the Debye approximations for the Bessel functions and the 
approximations (A 99) for the Legendre functions. We then get expressions of 
the type: 

; Oe 

Citin = 0 95 (1+ F) [Skee + Rkal. | 
(32) 


2 


») an a Oo » @) | 
CLo2 = @ | Si + Rxo2 — (1 + a) (Sxo2 + re) | | 


Here the S’s mean the contributions from terms with | < 9 — Vel and the 


R’s from the terms in the interval ae Vol =l>oe—kh 
The numerical computation of SY.,; »=1 or 2; meets no theoretical diffi- 


culties but is, of course, a very tedious piece of work. As a consequence of 
3 


the absence of saddle points in the interval (1; 9 — [Vo]) for 1 we must expect 
the sums Syoz to be of the same magnitude as a single term in the series. 
The calculations confirm this supposition; they have been carried out for o = 20; 
40; 60; n = 4/3; and a dozen values of y= 00. 

This is quite in accordance with the result obtained with the saddle point 
method. A fairly simple calculation shows that the asymptotic value thus afforded 
is zero of the contribution from the decisive point 1=0. The result also 
agrees with the geometrical optical picture given in A. A term with a certain 
value of J corresponds to rays with an angle of incidence a defined by: 


Sina : ce . (33) 


‘An upper limit for an exact calculation lies, for n=4/3, at about ep=18. Above this 


value the numerical work increases enormously, as we must then compute the Bessel func- 
tions ourselves. 
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We have therefore obtained the quite natural information that, for low values 


of a, the reflected part of the ray does not contribute to the intensity in the 
forward direction. 


There remains to calculate the influence from the few terms with 1 = 9 — Vol. 
Introducing the approximations (A 99) for the Legendre functions, we get: 


1) . A & 1)\0¢0 
Roz = ie F yoo: ; ly — (1 — 2) ? | ‘ 
r=1 Y | 
7 . mee (34) 
to = Froo(1—" : \Iofu— (7-9) a1;| 
y=1 ; 
n=o—l; q=[Vol. | 


It would be a great advantage to write these sums with the dependance on O 
outside the summation sign. Gegenbauer’s addition formula for Bessel func- 
tions! offers a possibility for this; for Jy and J, it runs: 


Jo (t) =Io(u) Jo(v) + 2 >) Im(u) Im (v) cos ma;| 
ee (35) 
Ns sea Ss Sete tat 2008 a | 
t ea ate v d cos a 
where: 
=u? —2uvcosa+ v*. (36) 


We observe that we can write the argument of the Bessel functions approxi- 
mately: 


27 — 4) (go 4+ 4)? 
- Introducing the new notations: 
cae 
pice! =a | 
; | (38) 
“a be 9 a 12 : 
erm = 3 (25 4)" (A) ro 
_we then get 
Fin = 2AM ,,, 4 g2OIW,,, | 
use (39) 


‘02 = [Jo (x) Jo(y) — 242 (2) Ja(y)] sx0 + 2 J, (x) Jy (y) sx1 + 
+ 4J5(x) J2(y) sxe; 


neglecting terms with a factor Jm(x) for m= 3. 


1 See Watson, Theory of Bessel Functions, p. 362. 
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We must now find appropriate expressions for Fyo2 as functions of n, 0 
and o. Leaving aside the case n ~ 1, we can still use (A 27) for qi1(n@), giving 


Vn? —1 y—+4 | 
a 4 : 40 
re n & e(n® — 1) A 
We further use the recurrence formula: 
; 1+1 i 
Ci(e) = ; Ci(o) — Cr41(@); (41) 


and thus get, omitting the argument 0 of the Bessel functions: 


en Let 4in 
HRY oleae Ty (42) 
Pua= Hits eee 
He, xo 2 
Here T, and T, denote: 
P= (nae hho Vy? 1 + ee ene 
; Ay Spe ; n@ Vn?—1 0 = 
es ; n “ . l P ( ) 
+8 — vena = 
To= ae — 1 + iV? -1 + ~ + ; 
; Ais 0 Vn? — Q | 


To compute the numerical values of the Hankel functions we use the Debye 
formula for the transitional region’ to get H®), (0): 


HH}; (0) = 0.246163 z F 0.062214 2? F 0.001026 z* + 0.000352 2 — 

— 0.000044 2? +... + 
+ i[— 0.426366 2 F 0.107758 22 + 0.001777 z4 + 0.000609 2° + { 
(ee 


and the recurrence formulae for cylinder functions to get the rest. 

Accordingly we calculate the quantities R via the sym’'s. The result is now, 
that in the forward direction the reflected radiation (provided the interpretation 
given in A is maintained) falls far short of the refracted. As an example we 
take y= 3.8; 0 =60. The moduli of both Co and Cyo2 are about 0.27 0. 
This is less than + of the refracted radiation in the same direction. 

This result does not agree with Mecke’s. His result is considerably larger. 
The difference is not surprising. The impossibility of using the saddle point 


* See Warson, Theory of Bessel Functions, p. 245. 
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method here must indicate that geometrical optics can describe the phenomenon 
but inadequately. 

There is another point worth noticing in the numerical results for Cos. The 
real part of Co, is much less than the imaginary; this corresponds to a do- 
minance of the real part of the field components and agrees with geometrical 
optics, 1.e. the expressions in Mecke’s paper yield the same phases. 

A grave objection can be raised to the calculations performed in this para- 
graph. It bears upon the carrying out of the summation only to 1=o—1. 
The next higher terms must have some influence on the phenomenon, in any 
case somewhat larger than the contribution from the reflected radiation. The 
author has persisted in performing the calculation to this fixed limit wholly 
from physical reasons. For the terms above this limit we lack a correspondence 
in the geometrical optical picture; to calculate the influence from them must 
therefore be performed from the original expressions, before the splitting up 
leading to the physical interpretation. Apart from this, the calculation can be 
carried out in the manner shown in this paragraph for the quantities R, but it 
is tedious. We will come back to this question in the next paragraph. 


4. On the positions of the minima 


We can actually use the expressions (4102) for the diffracted radiation 
unaltered. The only point to discuss is the passing from (A 100) to (A 102), 
but the exchange, if performed, is unimportant. We will in this paragraph be 
interested in the positions of the minima of intensity. Numerical checking shows 
that even for 9 = 25 the approximation (A 102) is safe. 

From (A 102) and (27) we can see that the sphere has no polarizing effect 
on the refracted and diffracted radiation in the forward direction. Thus we 
need consider one component of EF only. To get the simplest possible de- 
scription we will furthermore neglect the influence from the reflected radiation ; 
we can do this for two reasons, firstly the smallness of the terms Co., secondly 


ate ay 1 
the fact that the reflected radiation has a phase differing about a from that 


a 


of the diffracted radiation. 
Then we can write: 


E501 se E51 — i— cos gp -e thor (Foy ae Fs); (45) 
where: 
J1(y) 
Bee e : > | 
2 (46) 
. (n —1)(n + 1)? 


So we have to seek minima for: 


I=|Fo + Fil (47) 
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The first step to do this, the tabulation of y as a function of 6, reveals the 


important relation: 
l~l=1 (48) 


already utilized in the preceding paragraph. The difference for 6 < 1.5 never 
exceeds 2%, and the influence from the performed approximations must be 
considered larger. The physical meaning of (48) is that in the forward direc- 
tion, up to the angle giving the second minimum, the intensity of the refracted 
radiation is practically constant, but that the phase changes. 

On the other hand, the phase is a much more rapidly varying function of 
o, this being a consequence of the factor e~?'?”—), In comparison with the 
rapid fluctuation mentioned we can, for large particles (corresponding to low 
values of £), consider y independent of ©, using its value for 6 = 0: 


ke=0 = ils (49) 


Then we get the following simple picture. The minima are situated somewhere 
around the values given by (1a) in the intervals defined by: 


aye 


2 


(50) 


| 1 Qn 
ral 
Y 


i (ae nm + 1)? 


For increasing @ or 7 these intervals decrease; we also see that the observed 
minima are situated on one side of ym (as given by (1 a)) when m is odd, and 
on the other when m is even. Thus, in determining the radius of a large 
particle from the first two minima, using the »ideal> eq. (1 a), we get one 
value that is too large and another too small. 

It may interest to see the effect of an application of this rather rough 
approximation on an experimental material. We choose the measurements on 
mists and clouds given in the extensive treatise of KOHLER already cited. 
ery 
A 
ideal eq. (1 a) (7; = the value of the radius given in Kohler’s treatise) we try 
to find a better one, 0), being the least zero of the equation: 


J 3.8 aia 
ne ss2) ] 2n? 


Tt Sy eee ee ee 


Starting with the value 0, =~ obtained from the first minimum with the 


G1 3 g39 
01 


Afterwards we check our result by computing the effect of the corresponding 
correction on @g, obtained from the second minimum. From the table 4 in the 
cited treatise we get for the series 15—26!, that in 8 cases out of 9 the 
correction goes in the right direction; in three other cases the correction is 
small. An investigation of the series 69—90' in table 7 (measurements below 
0° C) gives the same result. Leaving aside 6 series in which the mean error 
exceeds 2%, we get a positive result in 11 cases and a negative in 2; in three 


' We only consider series with 7 > 11 u. 
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Fig. 1. 
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cases the change is small. This agreement is surprisingly good; in an investiga- 
tion on natural mists and clouds we must expect significant errors from un- 
foreseen sources, e.g. lacking homogeneity in the mist, or changes of its state 
luring the observations. Further, our theory has been based on the assumption 
of a monochromatic radiation. In his observations Kohler uses white light, in 
she calculations applying a mean of the wave lengths. Maintaining that the 
vgreement is better than could have been expected, we must at the same 
ime remember, that the comparison does not conclusively confirm our theory. 
[The uncertainty of the used material is too large. 
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Pig.c2: 
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We can find another conclusion drawn by Kéhler and supporting our simple 
picture. If we investigate a great many spheres, varying in size in an irregular 
manner and if we compute the values 7, and ry from the first two minima by 


ree , Neel & 
the ideal formula, then our description requires the mean of sae to be very 
2 


small. This is exactly what Kohler gets and he interprets it as an effect of random 
errors only. 

We can also compare with the experimental work of Mgcxe, loc. cit. On 
pp. 490 & 491 he gives a qualitative description wholly agreeing with what 
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has been said here. The spheres measured by Mecke are smaller than those 
discussed above. Moreover, it is to be regretted that he has only used optical 
methods to measure the radii of the liquid particles. Otherwise it would have 
been possible to make more than qualitative comparisons. 

When the radius is less than about 10 « we can no longer use (49) in (47) 
and the numerical work to get a substitute for it must be performed separately 
for each value of @ and ». The author has carried out these calculations for 
o=20 (5) 60; n=4/3; and the result is shown in diagrams 1 and 2. 
We see that the intensity curves for @ = 20 (10) 60 belong to a totally diffe- 
rent type than for 9@=25 (10) 55. This is a consequence of the fact that 
when o increases by 10, the angle 2e0(m— 1) increases by about 2 when n 
is 4/3, but when 0 goes from 20 to 25 the growth is only a little more than z. 
Then, the curves for o = 20, 30,... are so as to say homologous, the inter- 
ference between the refracted and the diffracted radiation being of about the 
same kind. For the other series, 90 = 25 etc., we have another set of homo- 
logous curves, and it is natural that the minima have gone to the other part 
of the interval given by (50). 

The following conclusions from the curves are worth mentioning. 


1. For increasmg values of @ the distance from the minima to the ideal 
ones decreases, if two homologous values of @ are considered. 

2. The first two minima lie on each side of the ideal minima in all the 
cases investigated. Hence, we may draw the conclusion that it is at least a 
very ordinary case, that the ideal equations (1a) give one value above and 
one beneath the correct value of the radius. 

3. It may happen, as for @ = 25 and 35 in the examples given, that the 
two minima run together to one broad minimum. When this occurs, the deter- 
mination of the size of the particles by measurements on the minima of 
intensity becomes impossible. Therefore it is not advisable to judge merely 
from the appearance of the phenomenon whether formula (1 a) can be applied 
or not. The diagram for 9 = 20 shows a case with »normal» intensity distribu- 
tion, but the deviation from the ideal curve is as large as for o@ = 25. 


Mecke also produces intensity curves of about the same kind: the omission 
of a factor » in F, does not, of course, change the essential character of the 
curve. But quantitatively the differences are large. The intensity in the 
secondary maxima are about doubled in the series @ = 20 (10) 60; and in the 
other series we get anomalous effects for larger drops than Mecke. Even the 
curve for 9 = 45 indicates that it must be difficult to observe two separate 
minima in that case. 

* 


We will now sum up what we can say as to the possibility of determining 
the radius of the sphere. For water (n = 4/3) it seems possible to obtain two 
values, one on each side of the correct one, whenever two separate minima are 
visible. This is always the case when @ exceeds 50. But this determination 
does not give a value as accurate as would allow the experimental data alone. 
To reach the best possible value we cannot restrict ourselves to observing the 
position of the minima only. It must be emphasized that the investigation 
peformed in this paper has used approximations that are perhaps very severe 
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for 0 <50. Especially this is the case for increasing 0, that is in the region 
of the second minimum. 

We will therefore take up the method suggested by Mecke, and using the 
colours of the phenomenon for @=:0 (the aureoles). The determination of 
these colours must be performed by means of Maxwell’s colour triangle in some 
modified form; the method used by Ként~ER* seems to be appropriate. As our 
expressions for the refracted radiation differs from Mecke’s, there must also 
result differences in the determined colours. Especially, we must expect coloured 
central maxima for larger drops than foreseen by Mecke. This is in accordance 
with the results of Kohler; he has obtained coloured aureoles already for 
# eolDbas 

In this paper we do not intend to carry out these calculations as they mean 
a great deal of numerical work. To be on the safe side, we can no more 
neglect the contribution from the terms Cy:, nor can we avoid a discussion of 
the terms with J surpassing 9. All necessary means for carrying out the rather 
tedious calculations are indicated in the preceding paragraph. It may be left 
to someone acquainted with the experimental side of the problem to judge, 
whether the importance of a solution justifies the necessary numerical pains. 


I wish to express my gratitude to Professor Ivak WALLER for suggesting 
this investigation and for his continued interest during the work. I am also 
indebted to Professor Hitpinc KOxHuLER for helpful discussions of the experi- 
mental side of this problem. 
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Errata in A 


P. 3, line 6, delete zero; p. 4, line —12, for Sinclair and Houghton, read 
Stratton and Houghton; p. 11, line 14, for sin & in the denominator read I sin 3; 
3 


| 1/2 % 
p. 20, line 2, . for / a read V/ m p. 21, line —9, for (@—%,), read 
2 


be 


3 er A 


‘ ; , 9 ; A : 270 2 ~ : - Vee 
(9—0)V; p. 21, Ime —7, fon are read 9 > PD 23, line —5, after Su insert 


FON (10) 9 (0) pyqq — FO (to) 9 (ta), BLE (aa) a) 

[P(e]? 4[ F(a)? ——-12-[F @o) FB 
p. 34, line —10, for 3, read O; p. 35, lines 26—27, for 0? read 0; delete but 
zero for 0~ x. 


(—1)!; p. 26, line 14, for — 
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